THE POTENTIAL {iz) m GENERATES REAL EIGENVALUES ONLY, 
UNDER SYMMETRIC RAPID DECAY CONDITIONS 



K. C. SHIN 

Abstract. We consider the eigenvalue problems — u" (z)±(iz) m u(z) = Xu(z), m > 3, under 
every rapid decay boundary condition that is symmetric with respect to the imaginary axis 
in the complex z-plane. We prove that the eigenvalues A are all positive real. 

o : 

^ \ Preprint. 

3 : 

■ 1. Introduction 

o : 

For integers m > 3 fixed and 1 < £ < m — 1, we are considering the non- Hermit ian 
^-i | eigenvalue problems 

1 '• Id 2 1 

m- (1) HMz):= —— + (-l) e (iz) m u(z) = \u(z), for some A G C, 

(N ' [ dz 2 

r- ; 

O with the boundary condition that 

O ' n £+1 

(2) u(z) — > exponentially, as z — > oo along the two rays argz = —— ± — ZjT^ 71 ' 

If a non-constant function u along with a complex number A solves flU) with the boundary 
condition (0), then we call u an eigenfunction and A an eigenvalue. It is known that any 
solution of (0) is entire (analytic on the complex plane). And it is also known that every 
solution of (0) is either decaying to zero or blowing up exponentially as z tends to infinity 
along any ray {zeC: arg z = const.}, except along m + 2 critical rays where the transition 
between decaying and blowing-up sectors might occur |T3"| , §7.4]. Along these m + 2 critical 
rays, any non-constant solution is decaying algebraically [ |T3"| , §7.4]. We will explain these 
asymptotic properties of the solution in Section 0. 

Before we state our main theorem, we first introduce some known facts about the eigen- 
values A of Hg, facts due to Sibuya [|19| and Hille ji~3"| . 



Proposition 1. For 1 < I < m — 1, the eigenvalues A of Hi are the zeros of a certain entire 
function of order | + — e (|, l) . In particular, the eigenvalues have the following properties. 

(i) Eigenvalues are discrete. 

(ii) All eigenvalues are simple. 

(iii) Infinitely many eigenvalues exist. 
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For our purposes, we will need to examine the proof of this proposition in some details. In 
Lemmas |7| and |8|, we prove that the eigenvalues are zeros of certain entire functions of order 
~ + — G f~, l). Then the claims (i) and (iii) are consequences of the Hadamard factorization 
theorem, while the claim (ii) is due to Hille |13|, §7.4]. 

In this paper, we will prove the following main theorem regarding the positivity of the 
eigenvalues. 

Theorem 2. The eigenvalues A of Hi for integers 1 < £ < m — 1 are all positive real. 

The eigenvalues of Hi are the same as those of H m _£, as we show in the proof by reflecting 
z i — > —z. The case £ = 1 of the theorem is due to Dorey, Dunning and Tateo f9|, and we 
use this in our proof. Also when m is even and I = one can see that Hi is Hermitian on 
the real line, and hence A 6 R. In all other cases Theorem ^| is new, and provides the first 
result of its kind for boundary conditions that neither cluster near the negative imaginary 
axis nor lie on the real axis. We will explain how Theorem ||] covers every symmetric rapid 
decay condition later when we discuss admissible boundary conditions in Section ||. 

For the rest of the Introduction, we will mention brief history and give some background 
about our main problem and our method of proof, and then in Section |2|, we will introduce 
work of Hille |13|, §7.4] and Sibuya |19| about some properties of solutions of (p. In Section 



[3], we establish an induction step on £, which is the key element in our proof of Theorem 0. 
More precisely, we will prove that for 1 < £ < — — 1, every eigenvalue A of H i+ i is positive 
real if all eigenvalues a of Hi lie in the sector | arg<x| < in the complex plane. In Section 
f|, we will outline a proof of the induction basis I = 1, that is, eigenvalues of Hi lie in the 
sector | arg A | < (in fact argA = 0; see |J and also [fTS[| ). We then prove Theorem 
|2] by induction on £ and the reflection z i— > —z. In Section |5|, we discuss some hopes and 
challanges in extending our method to more general polynomial potentials. Finally, in the 
last section, we mention some open problems. 

History and overview of the method. In this subsection, we introduce some earlier 
work related with our main result, Theorem |2|. Also, we discuss our method of the proof of 
Theorem |^. 

The Hamiltonians with the potential ±(iz) m have been studied in many physics and 



mathematics papers, either under a boundary condition on the real axis f| |20[, u(±oo+0i) 
0, or under the boundary condition (|2|) with £ = 1 |T|, p|, |18 . 



Simon |2(J and Caliceti et al. studied the Hamiltonians — -|C + x 2 + j3x m on the real 



line, where (3 G C — K._, m — 3,4,5, ... , and they proved compactness of the resolvent and 
discreteness of spectrum. Regarding the reality of eigenvalues, Caliceti et al. showed that 
eigenvalues for — ^ + x 2 + f3x 2n+l are real if (3 is small enough. 



Recently, a conjecture of Bessis and Zinn- Justin has been verified by Dorey et al. ||, (and 



extended by the author fl8|). That is, the eigenvalues A of 
d 2 



(3) 



- a(izy + (3z 2 



u(z) = Xu(z), u(±oo + Oi) 
verified for the case f3 



for «6l- {0}, peR, 
0, and later the author extended 



are all positive real. Dorey et al. 
for the case /?6R@. 

In fact, Dorey et al. have proved more. They studied the following eigenvalue problem 

d 2 



(4) 



dz< 



- (iz) 2M - a^zf 1 - 1 + 

z l 



u(z) 



\u(z), 



under the boundary condition (EJ) with £ = 1, and M, a, I being all real. They proved that for 



2l + l\, 
1 = 0), 



M > 1, a < M + 1 + |2/ + 1|, the eigenvalues are all real, and for M > 1, a < M + l 
they are all positive. A special case of (f§) is the potential iz 3 (when M = |, a - 
which is the (3 = version of the Bessis and Zinn- Justin conjecture, but their results do not 
cover the (3 ^ version. 



Later, the author [[18] studied the following eigenvalue problem. 
d 2 



dz 2 



+ [iz) 



+ aUizT' 1 + a 2 (iz) m - 2 + 



MIZ 



U(Z 



Xu(z) 



with the boundary condition (EJ) with t = 1, where G M. for all k. He proved that if for 
some 1 < j < y we have (j — k)ak > for all k, then the eigenvalues A are all positive real 
|18| , Theorem 2]. This covers the full BZJ conjecture. 

The proof of our main theorem, Theorem [|, has four parts. The first part follows closely 
the method of Dorey et al. ||, 0, developing functional equations for spectral determinants, 
expressing them in factored forms and estimating eigenvalues by Green's transform. The 
second part establishes an induction step on 1 < £ < y an d estimates eigenvalues by 
Green's transform again. In the third part, we use the result of Dorey et al. |10| that says 
the eigenvalues of Hi are all positive real (see also Theorem 2 in [[18]). This induction basis, 
along with the induction step established in the second part, proves our main theorem for 
1 < £ < y. Lastly we use the reflection z i— > — z to cover < £ < m — 1. Of course both 
this paper and p, |H], [L8[] are indebted to the work of Sibuya |19 . 



PT-symmetric oscillators. The above Hamiltonians are not Hermitian in general, and 
hence the reality of eigenvalues is not obviously guaranteed. However, these Hamiltonians 
share a common symmetry, the so-called "PT-symmetry. A PT-symmetric Hamiltonian is a 
Hamiltonian which is invariant under the combination of the parity operation V{: z \— > —z) 
(an upper bar denotes the complex conjugate) and the time reversal operation T(: i \— > —i). 
These PT-symmetric Hamiltonians have arisen in recent years in a number of physics papers, 
see P, 0, 12, pl| [16, ^l|] and other references mentioned above, which support that some 
PT-symmetric Hamiltonians have real eigenvalues only. The work of Dorey et al. || and 



the author |lq| , and the results in this paper, prove rigorously that some PT-symmetric 
Hamiltonians indeed have real eigenvalues only. 

We also know that if H = — -py + V(z) is PT-symmetric and V(z) is a polynomial, then 
V(z) = Q(iz) for some real polynomial Q, because V(— z) = V(z) and so KeV(z) is an even 
function and ImV(z) is an odd function. Certainly (|I|) is a PT-symmetric Hamiltonian. 
Moreover, if u(z) is an eigenfunction of Hi with the corresponding eigenvalue A G C, then 
so is u(—~z), with the corresponding eigenvalue A. 

The exceptional symmetric boundary condition. We will explain admissible boundary 
conditions in more detail in Section |2]. After that it will be clear that besides (Q) the 
only other exponentially decaying boundary condition that is symmetric with respect to the 
imaginary axis is that m is even and u(z) decays to zero as z — > oo along the both ends of 
the imaginary axis. 

But if u(z) is an eigenfunction of —u"(z) — (iz) m u(z) = Xu(z) satisfying this "exceptional" 
boundary condition, then v(z) := u(—iz) is an eigenfunction of the Hermitian equation 

-v"(z) + z m v(z) = -Xv(z), v{±oo + Oi) = 0. 

So A G R. 

The existence of the eigenvalues in this case is clear from Proposition |1] applied to the 
above equation of v. 

2. Properties of the solutions 

In this section we will introduce some definitions and known facts related with the equation 
(H|). One of our main tasks is to identify the eigenvalues as being zeros of certain entire 
functions, in Lemmas |7| and |8|. But first, we rotate the equation ([]]) as follows because some 
known facts, which are related to our argument throughout, are directly available for this 
rotated equation. 

Fix the integer m > 3. Let nbea solution of ([1]) and let v(z) = u(—iz). Then v solves 

-v"(z) + [(-lY +1 z m + A] v(z) = 0. 

And the boundary condition (fj) of u becomes that v(z) — > exponentially as z — ► oo along 
the two rays 

,t + l 

arg^ = ± ~^ n - 

m + 2 

Throughout this paper, we will use the complex number 



uj = exp 



2-rri 



m+2 



When I is even, it will be convenient to rotate once more, letting w(z) := v(u^z) so that 
w(z) solves 

-w"{z) + [z m + u\]w{z) = 0. 
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Figure 1. The Stokes sectors Sj of ([5]) with m = 3. The dashed rays are the 

7T I 37T 

5 ' 5 



critical rays: argz = ±?, 7r 



Hille [[13], §7.4] and Sibuya [0] have studied the general equation of the form 

-v"(z) + [z m + P(z) + X] v(z) = 0, 

where P(z) is a polynomial of degree less than m. 

We proceed now to summarize work of Hille [IBL §7.4], and expand on some work of Sibuya 



for P = 0, that is for 

(5) ~v"{z) + [z m + A] v{z) = 0, where m > 3. 

Results of Hille. It is known that every solution of ([5]) has simple asymptotic behavior 
near infinity §7.4]. We will explain this asymptotic behavior using the following. 

Definition. Consider the equation 

(6) g"{z) + [b m z m + b m ^z m ~ l + ... + b lZ + b ] g{z) = 0, 

where bk G C for < h < m with b m ^ 0. Let 

2jn - arg b m 

9j = — , for jeZ, 

m + 2 

where we choose —n < arg b m < n. For j G Z we call the open sectors 

(7) Sj = {zGC: % < arg z < 9 j+1 } 

the Stokes sectors of (||). Also we call the rays {arg z = 9j} the critical rays. 
In particular, the Stokes sectors of @ are 

(8) g J = (,GC: (2j - 1)7r <arg,< (2j + 1)7r l, for j £ Z. 
w J \ m + 2 6 m + 2 J J 

See Figure [T| 

Now we are ready to introduce some asymptotic behavior of solutions of (0). 
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Lemma 3 §7.4]). Let v be a non- constant solution of (|3p (with no boundary conditions 
imposed). Then the following hold. 

(i) In each Stokes sector Sj, the solution v is asymptotic to 



(9) (const.)z 4 exp 



2 TO+2 

-z 2 



as z — > oo in every closed subsector of Sj 



m + 2 

(ii) If v decays to zero in Sj, for some j G Z, £/ien mtzsi blow up in Sj-i and Sj+%. 
(However, it is possible for v to blow up in many adjacent Stokes sectors.) Moreover, 
the asymptotic expression (jQ) is valid with the same constant for the three consecutive 
Stokes sectors Sj-i U cl (Sj) U Sj+i, where cl (Sj) is the closure of Sj. 

(iii) For each Stokes sector Sj, there exists a solution of that decays in Sj, and there 
exists a solution of that blows up in Sj. And any solution of can be expressed 
as a linear combination of these two solutions. 

Remark. In Lemma |||, we state some asymptotic behavior of solutions of (§]). We mention 
that Hille fl3], §7.4] studied more general equations of the form (|6]). In the general case, the 
corresponding Stokes sectors are given by (]?]). With these Stokes sectors, Lemma ||| (ii) holds 
though the asymptotic expression @ becomes more complicated due to more complicated 
potentials. 

Also, one can see that Lemma |3] (iii) implies Proposition [I] (ii). More precisely, if eigen- 
values were not simple, then equation (JJ) would have two linearly independent solutions 
satisfying the boundary conditions in fl2|). Hence any solution of flU) could be expressed as a 
linear combination of these solutions. So there would be no solution of (P that blows up in 
either of the two Stokes sectors containing the rays in (0). This contradicts Lemma ^| (iii). 
Therefore, all eigenvalues are simple. 

From now on, we denote the Stokes sectors Sj as the Stokes sectors of (§). 

Admissible boundary conditions. Notice, in particular, that the asymptotic expression 
([]) implies that for each j, v(z) either decays to zero or blows up exponentially, as z ap- 
proaches infinity in closed subsectors of Sj. Also, the asymptotic expression (H) implies that 
if v(z) — > as z — ► oo along one ray in Sj, then v(z) — >• as z — > oo along every ray in Sj. 
Likewise, if v(z) —>■ oo as z —>■ oo along one ray in Sj, then v(z) —>■ oo as z —>■ oo along every 
ray in Sj. 

Let u be an eigenfunction of H^. Then the above observation shows that when I = 2n — 1 
is odd, the boundary condition (§) for Hi is equivalent to having v(z) = u(—iz) decaying to 
zero as z — > oo along rays in S- n and S n (note that the rays argz = ±^±^7r are center rays 
of S- n and S n ). Also, when £ = 2n is even, the boundary condition (fj) for He is equivalent 
to having w(z) = u(—i(u>^z)) decaying to zero as z — > oo along rays in S~ n -i and S n . 



As we saw above, one need not choose the two rays being symmetric, as in (0), so long as 
they stay in the Stokes sectors that are symmetric with respect to the imaginary axis. Next 
we explain why every solution of (|5|) decays to zero algebraically as z tends to infinity along 
the critical rays. To this end, we examine the asymptotic expression Certainly, one can 
check Re ^z 12 ^ J = for all z on the critical rays, and hence it is not difficult to see that 
every solution of (||) decays to zero algebraically as z tends to infinity along the critical rays. 
(Incidently, the Stokes sectors Sj are the sectors where Re (^z^~^j keeps a constant sign.) 

One might wonder why we do not consider the eigenvalue problem 



d 



2 



- (-i)^y 



u(z) = Xu(z), 



dz 2 

under the boundary condition ([|). (Here we have the opposite sign, compared to (|l|), in front 
of (—l) e .) In this case, under the rotation v(z) = u(—iz), and w(z) = v(cu^z) if necessary, we 
see that the two rays in (|2|) map to two of the critical rays of (|3p . So we have algebraic decay 
of the solution. Thus if we require the eigenfunction be decaying to zero 'exponentially', 
there are no eigenvalues, while if we require the eigenfunction be decaying to zero merely 
algebraically, then every complex number A is an eigenvalue. And hence we have no interest 
in this case. 

Now we are ready to explain how Theorem |2] covers every symmetric rapid decay boundary 
condition. When I is odd, the equation ([I]) becomes —u"(z) — (iz) m u(z) = Xu(z). In this 
case, the negative imaginary axis is the center of a Stokes sector and the critical rays are 
argz = — | ± for integers 1 < k < The two rays in (Q) are not critical rays 

and they are, in fact, centers of the Stokes sectors. When £ is increased by 2, the rays in (§) 
move, away from the negative imaginary axis, to the centers of adjacent Stokes sectors. So 
Theorem |2] covers all symmetric rapid decay boundary condition for the potential —(iz) m . 
Similarly, one can see that when I is even, Theorem |2] covers all symmetric rapid decay 
boundary condition for the potential (iz) m . 

So far in this subsection, we have discussed all possible symmetric decaying boundary 
conditions, except the "imaginary axis" boundary condition discussed at the end of the 
Introduction. Next we briefly mention non-symmetric decaying boundary conditions. Let 
us consider, as an example, 

d 2 



dz* 



- (iz) m 



u{z) = Xu(z) 



under the boundary condition that u(z) decays to zero exponentially as z tends to infinity 

f + ^andarg^-f-^. 



along the rays argz = — | + and arg^ = — | We set U\{z) := u{uj *z), and then 



ui(z) solves 



dz 2 



+ (iz 



, rn 



U\(z) = uj Xui(z) 



under the boundary condition (|2|) with £ = 2. Then by the results in this paper, Theorem [| 
one can see that u~ 1 X is positive real. Hence A is not real. In general, if we impose a decaying 
boundary condition along the two rays in some Stokes sectors that are not symmetric with 
respect to the imaginary axis, then the eigenvalues are not positive real. 

Finally, we mention that the integer £ in (|TJ) and (fj) is the same as the number of Stokes 
sectors between the two sectors where we impose the boundary condition (0). 

Results of Sibuya. Next we will introduce Sibuya's results, but first we define the order 
of an entire function g as 

log log M(r,g) 



order(g) = limsup 



logr 



where M(r,g) = ma 1 x{\g(re' ie )\ : < 9 < 2-rr} for r > 0. If for some positive real numbers 
a, Ci, C2, we have M(r,g) < ci exp^r ] for all large r, then the order of g is finite and less 
than or equal to a. In this paper, we choose arg z a = a arg z for — 7r < arg z < n and a G R. 

Now we are ready to introduce some existence results and asymptotic estimates of Sibuya 
jjTjl . The existence of an entire solution with a specified asymptotic representation for fixed 
A, is presented as well as an asymptotic expression of the value of the solution at z = as 
A tends to infinity. These results are in Theorems 6.1, 7.2, 19.1 and 20.1 of Sibuya's book 
The following is a special case of these theorems that is enough for our argument later. 



Proposition 4. The equation admits a solution f(z, A) with the following properties. 

(i) f(z,X) is an entire function of(z,X). 

(ii) f(z, A) and f r (z, A) = 4-f(z, A) admit the following asymptotic expressions. Let e > 0. 



Then 



z-^(l + 0(z- 1/2 ))exp 
-z^(l + 0(z- 1/2 ))exp 



2 m + 2 

-Z 2 



m + 2 



2 m + 2 

-Z 2 



,111 j 

(iv) 



(10) 
(11) 

(12) 



m + 2 

as z tends to infinity in the sector \ &rgz\ < — e, uniformly on each compact set of 
the complex X-plane. 

Properties (i) and (ii) uniquely determine the solution f(z, A) of ffl). 
For each fixed S > 0, f and f also admit the asymptotic expressions, 

f(0,X) = [l + o(l)]A- 1 / 4 exp 

/'(0,A) = -[l + o(l)]A 1 / 4 exp 
as X tends to infinity in the sector | argA| < 7r — 5, where 

k = J (y i + t m - Vt^^j dt. 

The entire functions X \— > /(0, A) and X \— > f'(0, A) have orders | + — . 



KA2- 



KX^ + ™ 



Proof. In Sibuya's book [[HJ, see Theorem 6.1 for a proof of (i) and (ii); Theorem 7.2 for a 
proof of (iii); and Theorem 19.1 for a proof of (iv). And (v) is a consequence of (iv) along 
with Theorem 20.1. Note that properties (i), (ii) and (iv) are summarized on pages 112-113 
of Sibuya's book. □ 

The next thing we want to introduce is the Stokes multiplier. Let f(z, A) be the function 
in Proposition |j. Note that f(z, A) decays to zero exponentially as z — > oo in So, and blows 
up in S-i U Si. Then one can see that the function 

(13) f k (z,\):=f(u- k z,u- mk \), 

which is obtained by rotating f(z,u~ mk X) in the ^-variable solves ([5]). It is clear that 
fo(z, A) = f(z, A), and that f k (z, A) decays in Sk and blows up in S k -iUS k+ i since f(z, u~ mk X) 
decays in Sq. Then since no non-constant solution decays in two consecutive Stokes sectors, 
fk and fk+i are linearly independent and hence any solution of (|5|) can be expressed as a 
linear combination of these two. Especially, for some coefficients Cj^(X) and Dj }k (X), 

(14) fj(z, A) = C j>k (X)f k (z, A) + D jik (X)f k+ i(z, A), j, k e Z. 

These Cj ik (X) and Dj k (X) are called the Stokes multipliers of fj with respect to f k and f k +i- 
We then see that 

(!5) C hk {X) = -f r-- and D jtk (X) ~ 



w k , k+ i(x) J>KV ' w k ,k+i{xy 

where Wj ik = fjf k — fjf k is the Wronskian of fj and f k . Since both fj and f k are solutions 
of the same linear equation ([5]), we know that the Wronskians are constant functions of z. 
Since f k and f k+ \ are linearly independent, W k , k +i ^ for all k G Z. In the next lemma, 
we will show that the Wronskian W ktk +i{X) is constant, which is needed in the proof of our 
main theorem. 

Lemma 5. For each fcgZ, the Wronskian W k ,k+i(X) = — 2u>~ I ^~ k ~ 1 , which is independent 
ofX. 

Proof. Since f k (z. A) = f(u~ k z,u~ mk X), we get 

f k+ i(z,X) = f{u~^z,u- m ^X) 
= f k {u- x z^- m X) 
= f k (u~ 1 z,uj 2 X), 

using u m+2 = 1. Also we see that 

f> +1 (z,X)=u- 1 f' k (u- 1 z,u 2 X). 

Thus 

(16) W j+ i M1 (X) = u- l W hk {uj 2 X). 
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Next we compute 

W-loCA) = Uf - fofU 

= f(a,z, u- 2 \)f'(z, A) - f(z, X)cof(coz, u- 2 X). 

Thus as z tends to infinity in So for which above asymptotics are valid, we have 

(uz) 



m + 2 
2 



m + 2 
m + 2' 



m + 2 
(JZ 2 



W_i i0 (A) = -M~T2T(i + 0(z-*)) exp 

+ z 1 lu(luz) * (I + <J(z 2)) exp 

= -2a; 4(1 + 0(2 2 )) since u; 2 = — 1. 
Finally we see that 

(17) W_i, (o,A) = -2a;-T J 

since is independent of 2. Thus (IB) and ( |TTD complete the proof. 



m + 2 
2 



m + 2 

z 2 



m + 2 



m+2 

z 2 



□ 



In the next lemma, we will show that the Wronskians Vt^^A) for 2 < k < m have orders 



+ 



which is needed in establishing the induction step in Theorem |TT] and proving the 
existence of eigenvalues under various boundary conditions. This lemma is due to Sibuya 
pl , but we give a full proof here. 



Lemma 6. For each 2 < k < m, the Wronskian Wk,o(X) has its order | + i £ (~, 1J. 

Proof. First we look at 

W M (A) = f k {z, \)f' Q (z, A) - f Q (z, X)f k {z, A) 

(18) = /(0, ^ 2fc A)f (0, A) - /(0, A)o;- fe /'(0, ^ 2fc A), 

since the Wronskian is independent of z and so we can take z — 0. We then know that 
the Wronskian Wfc^A) has order less than or equal to | + — because by Lemma [| (v) 
each term in the right hand side of ( |I8"D has order | + — . So in order to show Wk,o(X) 
has order equal to \ + ^, it suffices to show that there exist c\ > and C2 > such that 



|Wfc,o( A )l > ciexp 



C2IA 



m 

2 T m 



for all large |A| along some ray. 
Next we examine the right hand side of along the ray 

47T 



arg A = 9 



IT — 



m + 2 : 



by using the asymptotic expressions fllCf ) and (|TT|) . Notice 9 G (— 7r + 5, 7r — 5), for all 
< 5 < Recall that the expressions (|10|) and (|TT| ) are defined for arg A G (— ir + 5, re — 5). 
So in using (|ToD and ([□]) for f(0,u 2k X) and f'(0, u 2k X), we will choose 



(19) 



arg(^ 2fc A) =9* = 9 



Akn 
m + 2 



2tt, for 2 < k < 



m + 2 



11 



so we have 
(20) 

Thus we see 



mn „ „ ran „ , m + 2 
— 7T < < 9* < 9 < < 7T, for 2 < & < . 

m + 2 " _ _ , 



m 



f(0,u 2k X) = 

f'(0,u, 2k \) = 
Hence from this along with ([Tj|), we get 



[l + o(l)]|A| *e~ i ~*' exp 



K\X\ 



1 , J_ ■ m + 2 , 
2 T m g 2m 



[1 + o(l)]|A|*e*^ exp [i^Al^e^' 



W M (A) = -[l + o(l)]|A| 1 /Vl|A|^e-^ 



exp K\\\^ + ™ (g*G + ^) e + e i (3+£ 



+ [l + o(l)]|Ar 1/4 e-^- fc l 



exp tf|A| 



— e' * +w V -i 



[l + o(l)]exp K\\\* + ™ [ e l '(l+m, 



gK^ + m, 



-2iw-5[l + o(l)] exp ^|A|*+- ( gHl+w)" + g*U+w, 
where the last step is by 9 — 9* = 2n — • So when arg A = 7r — we nave 



Re 



Re 



cos 



m + 2, 
2m 



9 + cos 



m + 2 
2m 



ft, 



> 0, 



where the last step is by 

So for 2 < k < ^y*-, since K > 0, there exists C2 > such that 



|W M (A)| >exp 



c 2 



for all large | A | on the ray arg A = 7r 



47T 



m + 2 



Thus the order of W kfi (X) for 2 < jfc < ^ is \ + i 

Certainly the Wronskian Wfc^A) for 2 < k < is blowing up in some other directions 
as well. If one find 9 and satisfying (p0[), the above argument will show that Wk,o(X) is 
blowing up along arg A = 9. 



Since f k (z, A) = f(cu k z,u mfe A), we get 



fk+i(z, A) 



/(u;-( fc+1 ^,a;- m ( fc+1 U) 
/ fc ( W -^,a;- m A) 



,m+2 



using " = 1. 

Next for < k < m, one can choose 9 = —% + an d ^* by (HH)> an d then follow an 
argument similar to the above to conclude that the order of M4,o(A) is | + Or one uses 
an index change to reduce to the case already considered. That is, 



,o(A) 



fc-(m+2) 

u; m+2 ^^o, m+ 2- fe (u;- 2(m+2 - fc) A) by (16) 



—u 



m+2—kxir ( — 2(m+2— k) \\ 
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Then since 2 <m + 2- k < m ^ L , we know the order of W m+2 -k,o{^ 2{m+2 ~ k) ^) is \ + — , and 
hence so is the order of Wk,o{X) for < k < m. This completes the proof. □ 

Further results of Sibuya; Identifying the eigenvalues as zeros of certain entire 
functions. We can relate the zeros of C_ n .„_i(A) and D_„ jn (A) with the eigenvalues of H#. 
First, we study the case when £ is odd, as follows. 

Lemma 7. Let £ = 2n — 1 be odd, with 1 < £ < m — 1 . Then the following are equivalent: 

(i) A complex number A is an eigenvalue of Hi. 

(ii) X is a zero of the entire function C_ ni „_i(A). 
(hi) X is a zero of the entire function D_ ri) „(A). 

Moreover, the orders of the entire functions C_ ni „„i(A) and D_ nn (A) are ~ + — G (|, l) . 

Hence, the eigenvalues are discrete because they are zeros of a non-constant entire function. 
Note that the Stokes multipliers C_ nin _i(A) and D_ nn (A) are called spectral determinants 
or Evans functions, because their zeros are all eigenvalues of an eigenvalue problem. 

Proof. Suppose that A is an eigenvalue of H 2n _i with the corresponding eigenfunction u. We 
let v(z) = u(—iz), and then v solves 

-v"{z) + {z m + X)v(z) = 0, 

and decays in S- n US n . Since f- n (z, A) is another solution of the equation above that decays 
in S- n , we see that f- n (z, A) is a multiple of v. Similarly f n (z, A) is a multiple of v. Then 
since 



(21) f- n (z, A) = C^ n>n ^(X)f n ^(z, A) + D_ n>n ^(X)f n (z, A), 



and since f n -i( z , A) blows up in S n , we conclude that C_„ in _i(A) = 0. 

Conversely we suppose that C_„ )n _i(A) = for some A G C. Then from (ETj) we see that 
f- n (z, A) is a constant multiple of / n (z, A). Thus both are decaying in S_ n U S n , and hence 
u(z) := f- n (iz, A) is an eigenfunction of H 2n -\ with the corresponding eigenvalue A. 

Similarly, since 



f- n (z, A) = C_ n , n (A)/ n (z, A) + D_ riin (A)/ n+ i(z, A), 



we see that A is an eigenvalue of H 2n -i if and only if D_ n n (A) = 0. 

Finally, since Wk k+i(ty is a constant by Lemma ^| and since W^o(X) for 2 < k < m 
has order | + — by Lemma || we see from (0) that C_ n) „_i(A) and D_ n „(A) are of order 
| + — G (|, l) . This completes the proof. □ 

Second, we can relate the zeros of C_ n _i )fl _i(A) and D__„_i jn (A) with the eigenvalues of 
H#, when £ is even, as follows. 



Lemma 8. Let £ = 2n be even, with 1 < £ < m — 1. T/ien i/ie following are equivalent: 
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(i) A complex number X is an eigenvalue of Hi. 

(ii) X is a zero of the entire function A i— > C_ nin (u; A). 

(iii) X is a zero of the entire function X i— ► D_ n _i jn (o;A). 

Moreover, the orders of the entire functions X i— > C_ n)n (a>~ 1 A) and A i— ► D_ n _i )Ti ,(a;A) are 

2 ~ m C ^2> V - 

Note again that the Stokes multipliers C_ n)n (k> A) and D_ n _i ;n (u;A) are spectral deter- 
minants or Evans functions. 

Proof. Let u be an eigenfunction for H 2n with the corresponding eigenvalue A. Then v(z) : = 
u(— iz) solves 

v"{z) + z m v{z) = Xv{z). 
Next we let W\(z) := v(u~^z). Then W\(z) solves 

(22) -w'{(z) + [z m + or 1 A] w 1 {z) = 0. 

One can then check that the boundary condition for H 2n becomes that w\ (z) — > as z — ■> oo 
in U if> n+ i. Then like we did for Lemma 0, using 

aT 1 A) = C^ n {u)~ l X)f n {z,u)~ l X) + D_ n>n (w~ 1 A)/ n+ i(2, a; -1 A), 

one can show that A is an eigenvalue of H 2n if and only if C_ n „(u; _1 A) = 0. 
Similarly, we let w 2 (z) := v(u*z). Then w 2 (z) solves 

(23) -w' 2 \z) + [z m + uX]w 2 {z) = 0. 

One can then check that the boundary condition for H 2n becomes that w 2 (z) — > as z — > oo 
in S- n -i U S'n. Then using 

/_ n _i(z,u;A) = C^ n ^ n {uX)f n {z,uX) + D_ n _ 1>n (uX)f n+ i{z, uX), 

one can show that A is an eigenvalue of H 2n if and only if D_„_i jn (u;A) = 0. And this 
complete the proof. □ 

Next we want an infinite product representation of the entire function A i— > /(0, A). Recall 
that the integer m > 3 is fixed. 

Lemma 9. The functions f(0, A) and f'(0, A) have infinitely many zeros Ej < and Ej < 0, 
respectively. And they admit the following infinite product representations : 

(24) /(0, A) = A 1 f[ (l - ; or some A x e C - {0}. 

(25) /'(O, X)=A 2 f[ (l- for some A 2 e C - {0}. 
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Proof. We know that f(0, A) and /'(0, A) have orders | + — G (0, 1) by Proposition £| (v), and 
hence by the Hadamard factorization theorem (see, for example, Theorem 14.2.6 on page 
199 of |T3), we know /(0, A) and f'(0, A) have infinite product representations fl24|) and ( p5|) 
where /(0, = and /'(0, = for all j E N. So in order to complete the proof, we 
need to show E < if /(0, £) = 0, and E' < if /'(0, £') = 0. 

Suppose /(0, E*) = or /'(0, i?*) = 0. By the definition, we know f(z, £"*) solves 

-f"(z,E*)+z m f(z,E*) = -EJ(z,E*), 

and decays to zero exponentially in So- (Note that — + x m is Hermitian on the positive 
real axis.) In order to show E* < 0, we multiply this equation by f(z, E*) and integrate over 
the positive real axis to get 

/*oo roo poo 

- f"(x,E*)f(x,E*)dx+ / x m \f{x,E if )\ 2 dx = -E, \f(x,E*)\ 2 dx. 
Jo Jo Jo 

Next we integrate the first term by parts, and use /(0, E*) = or /'(0, E*) = 0. Then clearly 

the left hand-side of the resulting equation is positive, and hence we conclude E* < 0. This 

completes the proof. □ 

Next, we will prove a symmetry lemma, regarding f(z, A). 
Lemma 10. For each A G C, 



(26) f(z,X) = f(z,X) and f'(z, A) = f(z, A). 



Proof. This lemma for z = is contained in Lemma 8 in [18], and in fact that is all we will 



need in this paper. A proof of this lemma is essentially the same as that of 0, Lemma 8]. 



So we omit the proof here. □ 

3. The induction step 

In proving Theorem |^, we will use induction on I. The induction step will be provided by 
the following theorem. 

Theorem 11. If 1 < £ < y — 1 and all the eigenvalues of Hi lie in the sector \ arg-| < 
then every eigenvalue of Hf + i is positive real. 

Proof of Theorem [77] ( Case I: £ = 2n — 1 is odd, with 1 < n < Suppose that all the eigen- 
values (Tj of H 2n -i lie in the sector | arg<x,-| < -^h- That is, zeros o~j of the entire function 
cr i-^ D_ n n (a) lie in the sector | arg<r| < (see Lemma ^). Then we want to show that 
each eigenvalue A of H 2n is positive real. 

Suppose u(z) is an eigenfunction of H 2n with the eigenvalue A. That is, u(z) solves 

-u"(z) + (iz) m u(z) = Xu{z), 
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and decays to zero as z tends to infinity along the two rays argz = — | ± (2n + ■ ^ e 

then let v(z) = u(—iz), and so v(z) solves 

(27) + z m v(z) = Xv{z), 

and decays to zero as z tends to infinity along the two rays argz = ±(2n + 1)^X2 ■ 
Next we let w(z) = v(u~^z) so that w(z) solves 

(28) -w"{z) + [z m + u^X] w(z) = 0, 

and decays to zero as z tends to infinity along the two rays argz = (2n + 2)^2, —2n^^ 
that are the center rays of S n+ i and SL n . 
Then we examine 

or equivalently, 

(29) f(u n Zl u- 2n - l \) = C^ n (LO- 1 X)f(u- n z,u 2n - 1 X) + D^ n>n (u;- 1 X)f(cu- n - l z,^ n+1 X). 

We see that C- njn (uj~ l X) = by Lemma ||, and hence D_ njTl (a; _1 A) ^ 0. 

Next we will show that |D_ nj „(a; _1 A)| < f when A is an eigenvalue of H2 n with ImA > 0. 
To this end, we evaluate the equation (p9|) and its differentiated form at z = to get 

(30) fiO^-^X) = D_ n ,„(a;- 1 A)/(0 ) a; 2 " +1 A) and 

(31) w^W^A) = a;— 1 L>_ n , n (a;- 1 A)/ , (0,c l ; 2n+1 A). 

Then since A is an eigenvalue, so is A. Thus we have 

f^^~ 2n ~ l X) = D_ ri>n (a;- 1 A)/(0,a; 2n + 1 A) and 

^/'(o,^" 2 ™- 1 !) = ^-"-^^-^/'(oV^A). 

Then we take the complex conjugates of these, and apply (PU|) at z = to get 



/(0V" +1 A) = D_ n , n (u;- 1 A)/(0 ) a;- 2n - 1 A) and 



cu- n f(0,u 2n+1 X) = u n+l D^ n {^X)f{()^~ 2n - 1 X). 
So these along with (|3~0| ) and Q3l|) imply 



(32) D_ n , n (o;- i A)Z'- n)n (a;- 1 A) = 1. 

Clearly the order of the entire function a \— > D_ n) „(cr) is | + — G (0,1). So by the 
Hadamard factorization theorem we have 

00 / \ \ 

D- n ,n(X) = B Y[ i^l -J , for some B G C - {0}, 

where the Oj are the zeros of D_ njn (cr), and so | argOj| < for all j G N, by hypothesis. 
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Thus 



(33) 



-D-n,n(^ _1 A) 



n 

i=i 



n 



X 



CJCTj — A 



< 1, when ImA > 0, 



since Im (uhjj) > for all j G N. Hence this along with (^) implies 

\D_ nin (uj~ A) I < 1, when ImA > 0. 

Since the non-constant entire function f(z, u 2n+1 X) solves 

-f"(z,uo 2n+1 X) + [z m + uj 2n+1 X}f(z,u 2n+1 X) = 0, 

we know f(0, u 2n+l X) and f'(0,u 2n+1 X) cannot be zero at the same time. Otherwise, 
f{z,u 2n+1 X) = for all z G C. 

Suppose that f{0,u 2n+1 X) ^ 0. Then we get 

, ,— 2n— 1 \ 
oo 1 A 

'A) 



1 > LD-n 



n 



1 - 



(j2n+l A 



n 



2n+l 



- A 



u 2n+1 Ej - X 



> 1, when ImA > 0, 



where the last inequality holds since Im (u 2n+1 Ej) < if < arg(cj 2 " +1 ) < ir, by Lemma 
(this is where we use 1 < n < ~). So we have 



(34) 



ILL^^A)! 



when ImA > 0. 



Similarly, when f'(0, u 2n+1 X) ^ 0, we get 
1 > l£> 



-A)i=n 

i=i 



1 w- 2 "-^ 


oo 

=n 


u 2n+l E' 3 - X 




u 2n+l E' 3 - A 


1 3 







> 1, when ImA > 0, 



where the Ej are zeros of /'(0, So we again have (0). Hence ( |33| ) along with combining 
(0) and (H) gives 

(35) I I -U^l | 



n 



A 



Since |argcrj| < ^2 f° r an J ^5 we know Im (cucr,) > for all j G N. Moreover, 
Im (uj<7j) > for some j since both cr,- and ~a] are eigenvalues of the PT-symmetric oscillator 
i?2n-i- Therefore, from (35) we conclude A = A, and so A is real. 

We still need to show positivity of the eigenvalues. The function v(z) solves (p7[) and 



we know A G R. Also, one can check that v(z) solves the same equation. Then since the 
eigenvalues are simple, v(z) and v{z) must be linearly dependent, and hence v(z) = cv(z) 
for some c G C. Since \v(z)\ and \v(z)\ agree on the real line, we see that |c| = 1 and so 
\v(z)\ for all z G C. That is, \v(x + iy)\ is even in y. From this we have that 



viz 



(36) 







d_ 

dy 



\v(x + iy)\" 



-21m yu'(x)v(x)j , for all 



x G 
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Let g(r) = v(re l °) for 2n^^ < 9 < (2n + 2)^g. Then g(r) satisfies 

e~ 2i6 g"(r) + e mie r m g{r) = Xg{r). 



We then multiply this by g(r) and integrate over r e [0, oo) to get 

/•oo 

-2i8 / „///„\ 173 J„ i „mi0 



g"{r)g{r) dr + e 



r m |a(r)| 2 dr = A / \g(r)\ 2 dr. 



Next we integrate the first term by parts and multiply the resulting equation by e %e to get 

r m \g(r)\ 2 dr = Ae i 



-e- l ^'(0)^(0) 



|<?(r)| dr. 



|#'(r)| 2 dr + e (m+1)l9 

The we use e~ ie g'(0) = v(0) and take the imaginary parts to get, for all < 9 < 2r ^3^ , 

/*oo /*oo roo 

sin 6? / |g'(r)| 2 dr + sin[(m + 1)0] / r m \g{r) \ 2 dr = A sin 9 / |g(r)| 2 dr. 
Jo Jo Jo 



We evaluate this at 9 



2tvk 
m+2 



< 9 



(2ra+l)7r (2n+2)?r 
m+1 m+2 



(2n+i)7r ^ h ave A > 0. Note that since 1 < n < we see that 

m+1 — — 4 ' 



□ 



Proof of Theorem [77] f Case JJ: i = In is even, with 1 < n < 1JL ^- ). Suppose that all the eigen 

That is, zeros of the entire function 



values r,; of Ho n lie in the sector |are;r| < 



2tt 
m+2" 



r i — > D_ n _ l n (c<jT) lie in the sector | argr| < Then we want to show that each eigen- 

value A of H 2n +i is positive real. 

First, we examine D_ n _i n (r). From Lemma |] (iii), we know that the zeros of t h 
D_ n _ l n (u;r) are the eigenvalues Tj of H 2n , which lie in the sector | argr| < ^jr^, by hypoth- 
esis. So we have 

D_ n _ hn (uT) = B x Y[ I 1 - - ) , for some B x e C - {0}, 
where Im (uij) > for all j G N. Hence 

u^X 



D_^i, n (A) = BiJJ(l 



Thus 

(37) 



-n— l,n 



(A) 



f -n-l,n(A) 



n 



i - 



1 - 



n 

3=1 



UTj - A 



CJTj — A 



< 1, when ImA > 0, 



since Im (utj) > 0. 



Next, we suppose that u(z) is an eigenfunction of H 2n +i with the eigenvalue A. That is, 
u(z) solves 



-u"(z) + {-l) 2n+ \iz) m u{z) = Xu{z), 



and decays to zero as z tends to infinity along the two rays aigz = — | ± (2n + 2)^g. We 
then let u(z) = u(—iz), and so ^(2;) solves 

(38) -v"{z) + (^ m + X)v(z) = 0, 
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and decays to zero as z tends to infinity along the two rays argz = ±(2n + 2)^^ that are 
the center rays of S n+ \ and £_ n _i. 
Then we examine 

(39) f- n -i(z, A) = C^ n (X)f n (z, A) + D^ ljn (X)f n+1 (z, A). 

We see that C_ n _i )n (A) = by Lemma [7|, and hence D_„_i >n (A) ^ 0. So we have 



(40) f(uJ n+l z, oo- 2n - 2 X) = D-^Wfiu-"-^, oo 2n+2 X). 

Then we will show that |D_„_ l ra (A)| < f when A is an eigenvalue of H 2n +i with ImA > 0. 



Suppose /(0, 



UJ 



-2n-2 



A) 7^ 0. Since A is an eigenvalue, so is A. Thus we replace A by A in 



(fiOT), and then evaluate the resulting equation at z = to get 

/(0,o;- 2 "- 2 A) = D^ n (X)f(0,u 2n+2 X). 
Then we take the complex conjugate of this and apply ( ^6|) so that we have 

f(0,u; 2n+2 X) = D^jX)f(0,^ 2n - 2 X). 

Combining this with (|40|) at z = gives 

(41) D_ n _i, n (A)£)_ n _i in (A) = 1, 

since /(0, uj- 2n ~ 2 X) ^ 0. 

Similarly, when /'(0, u;~ 2n_2 A) ^ 0, we get (fH|) again. 

The equation ([II]) along with the inequality in (|37|) implies |D_ n _i )Tl (A)| < 1 when A is an 
eigenvalue of H 2n +i with Im A > 0. So we get from (|40|) at z = 

u,- 2 "- 2 A 



(42) 



1 > ID. 



n— l.n 



<A)i=n 

i=i 



1 



n 

j=i 



UJ 



2n+2 E, - X 



> 1, when ImA > 0, 



where the last inequality holds since Im (u + Ej) < if < arg(w ) < n (this is where 
we use 1 < n < ^j^-)- Then like in the proof for the case I odd, we have |D_ n _i in (A)| = 1, 
which is also obtained when /(0, u)~ 2n ~ 2 X) = (and hence f'(0, uj~ 2n ~ 2 X) ^ 0.) Therefore, 
we conclude A = A, and so A is real, like in the proof of the case t odd. 

We still need to show positivity of the eigenvalues A. Recall that AgK, and the function 
v(z) = u(—iz) solves (j38|). Let g(r) = v(re ld ). Then g(r) satisfies 



-2i6 Jit 



g"(r) + e m0 r m g(r) = Xg(r). 



We multiply this by g(r) and integrate over r G [0, oo) to get 

-2iB 



g"{r)g(r) dr + e 



mid 



r m \g(r)\ 2 dr = X / \g(r)\ 2 dr. 



Since v(z) decays exponentially to zero as z tends to infinity in S n +i U S- n -i, we have 
integratibility for (2n + 1)^ < 6 < ( 2n + 3 )^- 
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Note that since Aet, the equation ([36|) is valid for this case as well. We integrate the 
first term by parts and multiply the resulting equation by e 10 to get 



e-V (0)<7(0) - e 



-to 



\9 



o 



>(r)\ 2 dr + e {m+1)ie / r m \g{r)\ 2 dr = Xe ie \g{r)\ 2 dr. 



o 



o 



Then we use e g'(0) = i>(0) and take the imaginary parts with using (^) to get, for all 

(2n+l)7r n ^ (2n+3)7r 
m+2 ^ P ^ m+2 ' 

/*0O /"OO /*oo 

sin 6? / |(/(r)| 2 (ir + sin[(m + 1)#] / r m \g[r) \ 2 <ir = A sin 9 I \g[r)\ 2 dr. 
Jo Jo Jo 

We evaluate this at 9 = ( 2n+1 ) 7r to have A > 0. Note that since 1 < n < ^r^-, we see that 

m+l — — 4 ' 



(2w+1)tt ^ _ (2n+l)7r (2n+3)7: 
m+2 m+l m+2 



□ 



4. Proof of Theorem g] 

In proving Theorem |^, our induction basis is provided by the following lemma that is due 
to Dorey et al. |l(J (see also [18|, Theorem 2]). 



Lemma 12. The eigenvalues A of H\ are all positive real. 
Here we will give an outline of the proof. 

Proof. Let A be an eigenvalue of Hi with the corresponding eigenfunction u(z). Then we set 
viz) = u(—iz), and hence viz) solves 



-v(z) + (z m + X)v(z), 



and v(z) decays in S-i U Si. 
Then we consider 



(43) f-i(z, A) = C- lfl (\)f (z, A) + £>_ lj0 (A)/i(*, A). 

So we see that C_i i o(A) = 0. Moreover, we know that from fllED and (|i"6|) , and Lemma |^, 

W-i,o(A) 



l^-i,o(A)| 



VFo,i(A) 

Next, we use infinite product forms of either (|4"3"D when /_i(0, A) 7^ 0, or its differentiated 
form when fL±(0, A) 7^ at z = 0. Then like in the proof of Theorem [TTJ, one can show that 
A > 0, since the hypothesis in Theorem [TT] is needed for showing |£)_ 1)0 (A)| = 1 only. □ 

Now we are ready to prove Theorem 0. 

Proof of Theorem ^. For integers 1 < i < y, the theorem easily follows from induction on £ 



along with Theorem 11 and Lemma 12. 
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So we assume y < £ < m — 1. Suppose A is an eigenvalue of ifg with the correspond- 
ing eigenfunction U\(z). Then we let u(z) := U\(—z), and hence u(z) solves —u"(z) + 
(— lY(— iz) m u(z) = Xu(z). That is, u(z) solves 

-u"(z) + (-l) m -\iz) m u(z) = \u(z), 

since (— I) 21 = 1. 

Next, we examine the boundary condition. It is clear that since ui(z) decays along the 
two rays aigz = — ~ ± ^^vr, u(z) decays along the two rays 

7T £ + 1 n (m-£) + l 

argz = —ty 1 7r = 7r and 

s 2 m + 2 2 m + 2 

7T £+1 7T (m - I) + 1 

are z = 7r 7r = 1 tt, 

6 2 m + 2 2 m + 2 

which is the boundary condition with m — £ in the place of £. Hence u(z) is an eigen- 
function of H m _e with the corresponding eigenvalue A. Since ^ < £ < m — 1, we see that 
1 < m — £ < y- So by induction with help of Theorem [TT| and Lemma |T^, we conclude 



A > 0. This completes the proof. □ 



Remark. Our method in this paper closely follows the £ = 1 method of Dorey et al. [10] and 



the author |Tj|. One big difference is that £ = 1 implies |D_i ) o(o')| = 1 for all a & C, while for 
1 < £ < m — 1, the corresponding functions a i— > D_„ n (u; _1 cr) in ( p9|) and a i— > D_ n _ l n (cr) 
in ( |39"D are entire functions of order ^ + ■ However, when A is an eigenvalue of Hg + i, under 
some hypothesis on the eigenvalues a of He we are able to show that |D_ ni „(a; _1 A)| = 1 for 
£ = 2n — 1 odd, and |D_ n _i ]n (A)| = 1 for £ = 2n even. This is the new and main idea in 



proving the induction step, Theorem 11 



5. Hopes and challanges in extending to more general potentials 

In this section, we mention some hopes and challenges in extending our induction methods 
to more general polynomial potentials. As mentioned at the end of the previous section, 
we used the induction on £ to show that |D_ nn (u; _1 A)| = 1 for £ = 2n — 1 odd, and 
|-D-n-i,n(A)| = 1 for £ = 2n even. Of course, if one can find another way of proving 
these identities it will help us find some new reality results of eigenvalues for more general 
polynomial potentials. 

In [1181, the author studied 

(44) Hfu(z) = 

under the boundary condition (|) with I — 1, where P(z) = axz m ~ x + a 2 z m ~ 2 + ■ — V a m -\Z 
is a real polynomial. I proved that if for some 1 < j < f we have (j — k)ak > for all k, 
then the eigenvalues A are all positive real. 



dz 2 



+ (-lY(iz) m + P{iz) 



u(z) = Xu(z), 
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One might wonder whether or not it is possible for 1 < i < m — 1 to combine our method 



in this paper and the method in author's earlier paper fllSfl , in order to study the reality 
of the eigenvalues of the ("PT-symmetric) equation ( [44]) under the boundary condition (|2|). 
The work of Hille [|IJ| §7.4] mentioned in Section [2] certainly covers the equation of the form 
([4*4]), and Sibuya [19| studied the more general equation ([44]) in a rotated form (even though 
in this paper we state only a special case, Proposition f|, of his results for simplicity). So one 
can check that the material in Section [2] has a natural generalization (see Section 2 in JIB 
for some generalization). 

Moreover, the spectral determinants in these cases, corresponding to those in Lemmas [1] 
and [|, are entire functions of the coefficient vector a : = (a m _i, a m -2, • • • , ai) £ C m_1 . Hence 
the eigenvalues Xj(a) are continuous in the coefficient vector a £ C m_1 . With a = we 
know the Xj(0) are all positive real by Theorem |2| in this paper, and so for each j £ N 
there exists a neighborhood of the origin in C m_1 such that if the coefficient vector a lies 
in this neighborhood, we have | argAj(a)| < j^rg. So if the coefficient vector a lies in the 
intersection of all such neighborhoods of the origin in C m -\ so that |argAj(a)| < ^ for 
all j £ N, then using the ideas in this paper one would prove, with a little effort, the reality 
of the eigenvalues, provided a £ IR m_1 with some sign restrictions on like for the £ = 1 
case above. Here we restrict a £ IR m_1 because this ensures that the eigenvalues appear in 
complex conjugate pairs or else are real, which is needed in the proof of the analogue of 
Theorem [11]. This along with inequalities | argAj(a)| < for all j E N proves identities 
that correspond to |D_ n) „(u; _1 A)| = 1 for £ = 2n — 1 odd, and |D_ ra _i jn (A)| = 1 for £ = 2n 
even. Once we have these identities, the remaining portion of the proof of reality of the 
eigenvalues follows closely the method in [[TBJ]. However, the author does not know whether 
or not the intersection of all such neighborhoods of the origin in C m_1 contains more than 
just a = 0. 

If the coefficient vector a £ M m ~ 1 lies outside the intersection mentioned in the previous 
paragraph, it seems that one should try to separately use induction on I. Unlike for the case 



P = 0, when P ^ 0, if we try to establish an induction step similar to Theorem [IT], the 
equation corresponding to ( "28]) in the proof Theorem [IT] for £ = 2n — 1 odd becomes 



-w"{z) + [z m - LO^Piu-h) + u^A) w{z) = 0, 



with w decaying in S n+ i U S^ n . So following the method in the proof of Theorem |IT|, we need 



to have | arg cr^- 1 < if o~j is an eigenvalue of H^Zi where P2 n ~i{z) := — u~ l P{u~^z) 
that is non-real. There is no reason to believe the eigenvalues of H^Lx (which is not VT- 
symmetric) are all positive real. So there could perhaps be some non-real eigenvalues. A 
challange comes from not knowing a method of estimating the arguments of the non-real 
eigenvalues. The infinite product expressions like (|33|) are very effective for proving the 



22 



reality of the eigenvalues in this paper, but are they still useful in estimating the arguments 
of the non-real eigenvalues? Perhaps not. 



In |I7| , the author studied the equation ( f44|) with m odd, under the decaying boundary 



condition at the both ends of the real axis, and proved the eigenvalues A lie in the sector 
largAI < ~r"o under some restriction on the coefficient vector, but still with some large 

I o i — m+2 ' ° 

a G lR m_1 . So we have a promising method of estimating the arguments of the non-real 
eigenvalues. But in the current case, we need to have a reverse induction on 1 < £ < y. One 
can check with a little effort that even for P = 0, this attempted reverse induction does not 
give us the desired identities | D_ n>n [u~ x A) | = 1 for £ = In — 1 odd, and |£L n _i in (A)| = 1 for 



£ = 2n even, that are needed in our proof of Theorem |TT| So one needs to show identities 
like these perhaps in a different way. 

Finally, a simple way of seeing why this reverse induction on 1 < £ < y fails is as 
follows. In establishing the induction step in Theorem [n], we restrict £ on the interval 
1 < £ < y — 1 because we have some technical difficulty that one can see from sentences 



right before equation fl34|), and right after equation (f42|). From this observation, it is not 
difficult to see that trying to establish a reverse induction step on 1 < £ < y is essentially the 
same as trying to establish the forwards induction step on | < f < m - 1 that our method 
is unable to provide. (In proving Theorem |^ for y < £ < m — 1, we used the reflection 
z I— > —z.) 

6. Conclusions 

In this paper, we proved that for each integer 1 < £ < m — 1, eigenvalues A G C of 

-u"{z) + (-l)\iz) m u(z) = Xu(z), 

under the boundary condition that u(z) — > exponentially, as z — > oo along the two rays 
argz = — | ± are all positive real. And we studied other boundary conditions. Due 
to asymptotic behavior of the solution u(z) near infinity, if the Stokes sectors that contain 
the two ray where we impose the boundary condition are symmetric with respect to the 
imaginary axis, then the eigenvalues are all positive real, except the case when the two 
Stokes sectors contain the imaginary axis. For all other boundary conditions, either there is 
no eigenvalue or the eigenvalues are not real. 

It will be interesting to consider the eigenvalue problem with more general polynomial 
potentials (— 1) (iz) m + P{iz) where P(z) is a real polynomial with degree less than m, 
under the boundary condition (|2D. 

It is known that some PT-symmetric oscillators with some cubic and quartic polynomial 
potentials have non-real eigenvalues [0, || [11], |I2|| . One would like to classify when VT- 
symmetric oscillators with polynomial potentials have non-real eigenvalues. 



23 



Acknowledgments. The author thanks Richard S. Laugesen for continued encouragement 
throughout the work. 

References 

[1] C. M. Bender and S. Bocttcher. Real spectra in non-Hcrmitian Hamiltonians having PT-symmetry. 

Phys. Rev. Lett, 80:5243-5246, 1998. 
[2] C. M. Bender, S. Boettcher and V. M. Savage. Conjecture on the interlacing of zeros in complex Sturm- 

Liouvillc problems. J. Math. Phys., 41:6381-6387, 2000. 
[3] C. M. Bender and A. Turbiner. Analytic continuation of eigenvalue problems. Phys. Lett. A, 173:442- 

446, 1993. 

[4] C. M. Bender and E. J. Weniger. Numerical evidence that the perturbation expansion for a non- 
Hermitian PT-symmetric Hamiltonian is Stieltjes. J. Math. Phys., 42:2167-2183, 2001. 

[5] C. Bernard and V. M. Savage. Numerical simulations of "PT-symmetric quantum field theories. Phys. 
Rev. D, 64,085010:1-11, 2001. 

[6] E. Caliceti, S. Graffi and M. Maioli. Perturbation theory of odd anharmonic oscillators. Comm. Math. 
Phys., 75:51-66, 1980. 

[7] E. Dclabaere and F. Pham. Eigenvalues of complex Hamiltonians with PT-symmetry I, II. Phys. Lett. 
A, 250:25-32, 1998. 

[8] E. Delabaere and D. T. Trinh. Spectral analysis of the complex cubic oscillator. J. Phys. A: Math. Gen., 
33:8771-8796, 2000. 

[9] P. Dorey, C. Dunning and R. Tateo. Spectral equivalences, Bethe ansatz equations, and reality properties 

in PT-symmetric quantum mechanics. J. Phys. A: Math. Gen, 34:5679-5704, 2001. 
[10] P. Dorey and R. Tateo. On the relation between Stokes multipliers and T-Q systems of conformal field 

theory. Nucl. Phys. B, 563:573-602, 1999. 
[11] C. R. Handy. Generating converging bounds to the (complex) discrete states of the P 2 + iX 3 + iaX 

Hamiltonian. J. Phys. A: Math. Gen., 34:5065-5081, 2001. 
[12] C. R. Handy, D. Khan, Xiao-Qian Wang and C. J. Tymczak. Multiscale reference function analysis of 

the PT symmetry breaking solutions for the P 2 + iX 3 + iaX Hamiltonian. J. Phys. A: Math. Gen., 

34:5593-5602, 2001. 

[13] E. Hille. Lectures on Ordinary Differential Equations. Addison- Wesley, Reading, Massachusetts, 1969. 

[14] E. Hille. Analytic function theory, Volume II. Chelsea Publishing Company, New York, 1987. 

[15] G. A. Mezincescu. Some properties of eigenvalues and eigenfunctions of the cubic oscillator with imag- 
inary coupling constant. J. Phys. A: Math. Gen., 33:4911-4916, 2000. 

[16] A. Mostafazadeh. Pseudo-Hcrmiticity versus PT Symmetry: The necessary condition for the reality of 
the spectrum of a non-Hcrmitian Hamiltonian. J. Math. Phys., 43:205-214, 2002. 

[17] K. C. Shin. On the cigcnproblems of PT-symmetric Oscillators. J. Math. Phys., 42:2513-2530, 2001. 

[18] K. C. Shin. On the reality of the eigenvalues for a class of PT-symmetric Oscillators. To appear in 
Comm. Math. Phys., 23 pages, 2002. 

[19] Y. Sibuya. Global theory of a second order linear ordinary differential equation with a polynomial coef- 
ficient. North-Holland Publishing Company, Amsterdam- Oxford, 1975. 

[20] B. Simon. Coupling Constant Analyticity for the Anharmonic Oscillator. Ann. Phys., 58:76-136, 1970. 

[21] M. Znojil. Spiked and PT-symmetrized decadic potentials supporting elementary 7V-plets of bound 
states. J. Phys. A: Math. Gen., 33:4911-4916, 2000. 

EMAIL CONTACT: kcshin@math.uiuc.edu 

Department of Mathematics, University of Illinois, Urbana, IL 61801 



